Magnon Modes for Thin Circular Vortex State Magnetic Dot 
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The magnetization in a magnetic microdot made from soft magnetic materials can liave a vortex- 
like ground state structure resulting from competition between the exchange and dipolar interactions. 
Normal mode magnon frequencies for such dots are calculated taking into account both exchange 
and magnetostatic effects. The presence of a low-lying mode as well as doublet structure with small 
splitting is demonstrated. Estimates of the mode frequencies for permalloy dots are obtained, and 
the possibility of experimental detection of such modes is discussed. 

PACS numbers: 75.60. Jp, 75.10.Hk, 74.40. Cx 



Nanostructures fabricated from soft magnetic materi- 
als are of interest both in basic and applied magnetics 
with potential applications including high-density mag- 
netic storage using an array of dot structures. BeI In the 
following we will consider the dynamics of thin mag- 
netic dots, with the important parameter in relation to 
the dot radius, R and the dot th ickness,!/ being the ex- 
change length, Iq = \/ Aj AttM'^ , where A is the inho- 
mogeneous exchange constant and JMs is the saturation 
magnetization. It has been showno that under certain 
conditions a vortex structure will be stable because of 
competition between the exchange and dipole interac- 
tions, and it is expected that these non-uniform states 
will drastically change the dynamic and static proper- 
ties of a dot. For this reason dynamic properties for thin 
{L <§; R) magnetic dots in the vortex state are studied by 
determination of the spin wave normal mode frequencies. 
The method used is based on the deteapiination of the 
of vortex-magnon scattering amplitude,lj where the non- 
local dipolar interaction is included by calculation of the 
magnetostatic energy arising from spin waves on the vor- 
tex background. A calculation of the discrete spin-wave 
spectrum leads to modes labeled by an azimuthal num- 
ber, TO, and a principal number, n. The magnon modes 
for the dot are compared with modes in the uniformly 
magnetized dot.B For vortex state dots a quasi-Goldstone 
mode (to = 1 , n = 0) with an anormalously small fre- 
quency appears and the doublets with to = ±|to| are 
split. 

For consideration of the vortex state it is natural to 
write the magnetization, M by use of angular variables in 
the polar coordinate system AI^ = Ms sm9 cos (p, My — 
Ms sin 9 sin ip, M^ = Ms cos 9. For thin enough dots with 
L smaller that the vortex core size the magnetization is 
uniform along the dot axis z, and the vortex state is given 
by the ansatz 



9^9ir),ip^X±Vo- (1) 
where r and x a-re polar coordinates in the dot plane. 



This distribution is typical for magnetic vortices in 
two-dimensional easy plane (EP) ferromagnet (FM), and 
the function 9{r) is determined Iw numerical solution of 
an ordinary differential equation!] This gives the results 
that in the center of the dot (r = 0) sin 6'(r) = and for 
r :s> Iq, the value 9{r) tends to 7r/2 exponentially. For 
the magnetic dot made of soft magnetic materials the 
crystallographic EP anisotropy is negligible, and the de- 
magnetization field, Hm is dominant. Assuming that the 
magnetization along the dot axis z is uniform, the energy 
is as an integral over the dot plane. 
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The sources of magnetostatic field H„i = — V$ are 
both volume and surface " magnetic charges" arising from 
divM, and from the discontinuity of the normal com- 
ponent of M on the surface. For determination of the 
demagnetizing field the vortex distribution (|^) is very 
simple because divM = and Mr\r=R for ipo = 7r/2, 
thereby giving neither dot edge nor volume contribu- 
tions and the sole source of field, for the ground state 
is Mz- For a thin enough dot this gives the contribution 
to the energy —MHm/2 = 27rAf^ which is an effective 
EP anisotropy forming the out-of-plane vortex structure. 
The magnetostatic potential from the edge and volume 
are determined by 
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'-Rdxdz (3) 



r=R 



where (f— r ~ r^-f r'^ — 2rr' cos(x — x') + (z — z')^. The 
magnon modes are investigated using the usual procedure 
of expansion in small deviations from the vortex state 
using 9 — 9o{r) + 'd and (p — x + 7r/2-|-^/;. It is convenient 
to define the new small variable, fi = — ■(/'sin^o so that 
the corrections to magnetization can be expressed in the 
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simple form, SM/Ms — "d ^xcos6'o — zsuiOq^ + ^if. The 
linear equations of motion for /i, -d can be written as 



where $ = $t, + '^edge , from (^), and G is the same d 
mensionless operator as for the vortex in EP magnets. 
Thus, in contrast to the case of EP magnets with lo- 
cal interactions^ the equations for and i) and /i become 
integro-differential. 

To proceed with the solution of (Q) the ansatz z9 = 
fmii^) cos(x'7i + cot), ji = gmif) sin(x»Ti -I- tut) is used. It 
is easy to see that for this ansatz the values of VM and 
M'r\r=R , as well as the magnetostatic potentials and 
^edge , are proportional to sin(x'Ti -I- cut). Then, one ob- 
tains the set of ordinary integro - differential equation for 
functions/m (r) , gm (r) instead of with partial deriva- 
tives. 

The edge contribution to Hm can be estimated from 
the integral (||), which diverges as L — > giving a non- 
analytic contribution to the system energy for small L. 
Using asymptotic techniques, the edge contribution to 
the magnetostatic field can be written in the form Hr = 
-fi{R,x)MMr), where F™(r) = (L/i?)(r/i?) as 
r —> 0, and reaches its maximal value independent of the 
small parameter L/ R, with Fm{r) —>■ 2TT , as r —>■ R near 
the edge. The maximal value of the volume contribution 
to Hr is small (the order of L/R) over all of the dot. 
Thus, the main contribution to the magnetostatic energy 
is determined by the field Hr in the region near the edge 
R — r < L. The leading (logarithmic) contribution to the 
energy in the parameter L/R for any m is determined 
from the value of g on the edge only, 



2TrL'^RImM^g^{R) In 
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(5) 



where Im — 2 for m — and Im — 1 for the other modes. 
Therefore, to lowest order in L/R , the contribution from 
non-local magnetostatic potentials can be accounted by 
use of the effective boundary condition for (?(r). 
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To find the frequencies of magnon modes in the low- 
est (logarithmic) order in L/R , one can solve (^) with 
$ = 0, find fi{r) and apply the boundary condition 
(^. The solution far from the vortex core, r ^ Iq is 
gm{r) = Jm{kr) + <7m{klo)Ym{kr) , where J™(z),y,„(z) 
are Bessel and Neumann functions, respectively, fc is a 
wave number, and am{klo) is the so-called scattering 
amplitude^ For the case of interest kR < 1 the value of 
klo £ Iq/R ^ 1 , and long- wave approximation is valid. 
In this case, the spin wave frequency uj = ujolok, where 



ujQ = 'iTT"fMs, and the asymptotics for (T,„(fcZo) found inE 
can be used. For estimates, we will use the values for 
permalloy, A^^^/M^ — 17 nm, that gives Iq — 4.8 nm 
and the characteristic demagnetizing field AttMs = 1.012 
T, leading to the following results. 

In the first approximation on Hq <^ 1 the scattering 
amplitude is small and gmi^) — Jm{kr). The frequencies 
of modes with m = \m\ and m = —\m\ are equal and 
are determined by ujn,m = ^n,\m\ = '^ir^Ms {nn^mh/R), 
where Kn,m is the n-th solution of kJ^(k) -f AJm(K) = 0. 
The value of K,n,m falls between the n -th root of the 
Bessel function im,n and the rt-th root of its derivative 
Jm,n> ~^ jm,n at A ^ m and k — > „ in the opposite 
case. Taking intaaccount <Jm, the value of which depends 
on the sign of to,EI produces the splitting of these doublets 
for the modes with m ^ 0. The scattering amplitude is 
maximal (linear in Jc) for the modes with m ± 1, with 
(Jm = 7rA:Zow/4 |m|,Q and the splitting of this doublet is 
the order of the next power on the small parameter Iq/R, 
uJm=-i — ^m=i ~ 6.16Li;o(^o/-R)^- Thus, for a typical R 
like 100 nm, this splitting is of order (0.5 2)GHz and is 
much smaller than the mean frequency, see Fig.l. 

Large values of the scattering amplitude leads to the 
most interesting feature of the vortex state dot frequen- 
cies, namely, the appearance of a low lying mode with 
TO = 1 and n = corresponding to vortex displace- 
ment from the center of the dot, the so-called trans- 
lational Goldstone mode CTGM), for which Ji{kR) ~ 
{klo)Ym{kR) and kR ^1.0 For vortex state magnetic 
dots the frequency of the TOM mode, shown in Fig. 2 , de- 
pends non-monotonically on R, ujtgm = ^o{lo / R)'^ {-^ ~ 
1)/(A-|- 1). Comparison of Figs. 1 and 2 shows that it 
has the same order of magnitude as the doublet split- 
ting for higher modes with \m\ = 1, and it much smaller 
than 1371 or smaller than frequencies for homogeneously 
magnetized dots. The frequency of TGM found from the 
effective equation of motion for the coordinate of the vor- 
tex X with some model assumptions about static vorte, 
energy shows different dependence on the dot radius,Q'l 
namely, u oc l/M^ in Ref . ^ and a; c>c 1/i? in Ref . ^. In 
contrast with,acl we do not need in any model assumption 
and our approach is based on the analysis of Eq. (^) . 

The modes with to = are singlets having cylindrical 
symmetry, and their frequencies for ti = and n — 1 are 
shown in Fig. 1. The lower mode with m = 0, n = 
can be considered as coupled oscillations of the vortex 
core size and ipo defined in (^. It can be observed in 
resonance experiments with an ac field perpendicular to 
the dot plane. 

Modes with \m\ > 1 have small values of am oc (klo)^, 
77 > 4 for different to with small doublet splitting, 
{ujn^-m — ^n,m) /'^n,\m\ ~ {Iq / R^ , that cannot be seen 
on the scale used for Fig. 1. The mean frequencies -ajn i^i 
of these modes for the vortex state dots have the same 
order of magnitude as for the homogeneous states modes 
calculated in.B 

In summary, it is shown that the dynamics of the mag- 
netic dot in the vortex state differs significantly from 
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the case of uniformly magnetized dots. The magnon 
mode corresponding to oscillations of the vortex position 
(TGM, m — l,n — 0) have the lowest frequency. The 
next mode (m = 0, n = 0) describes the oscillations of 
the vortex core size. These two modes can potentially be 
detected in resonance experiments, with ac - field parallel 
and perpendicular to the dots plane, respectively. Higher 
modes with to ^ forms doublets with maximal splitting 
less that 1 GHz. 
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FIG. 1: Mode frequencies as a function of dot radius R. 



FIG. 2: Frequency of the TGM mode versus dot radius R 
for different values of dot thickness L. 



